home *** CD-ROM | disk | FTP | other *** search
/ Cream of the Crop 26 / Cream of the Crop 26.iso / os2 / octa209b.zip / octave-2.09 / DLFCN.ZIP / dlfcn / examples / oregonator.cc < prev    next >
C/C++ Source or Header  |  1997-03-07  |  368b  |  18 lines
  1. #include <octave/oct.h>
  2.  
  3. #include <iostream.h>
  4.  
  5. DEFUN_DLD (oregonator, args, ,
  6.   "The `oregonator'.")
  7. {
  8.   ColumnVector dx (3);
  9.  
  10.   ColumnVector x = args(0).vector_value ();
  11.  
  12.   dx(0) = 77.27 * (x(1) - x(0)*x(1) + x(0) - 8.375e-06*pow (x(0), 2));
  13.   dx(1) = (x(2) - x(0)*x(1) - x(1)) / 77.27;
  14.   dx(2) = 0.161*(x(0) - x(2));
  15.  
  16.   return octave_value (dx);
  17. }
  18.